The High Redshift Integrated Sachs- Wolfe Effect 
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In this paper we rely on the quasar (QSO) catalog of the Sloan Digital Sky Survey Data Release Six 
(SDSS DR6) of about one million photometrically selected QSOs to compute the Integrated Sachs- 
Wolfe (ISW) effect at high redshift, aiming at constraining the behavior of the expansion rate and 
thus the behaviour of dark energy at those epochs. This unique sample significantly extends previous 
catalogs to higher redshifts while retaining high efficiency in the selection algorithm. We compute 
the auto-correlation function (ACF) of QSO number density from which we extract the bias and the 
stellar contamination. We then calculate the cross-correlation function (CCF) between QSO number 
density and Cosmic Microwave Background (CMB) temperature fluctuations in different subsamples: 
at high z > 1.5 and low z < 1.5 redshifts and for two different choices of QSO in a conservative 
and in a more speculative analysis. We find an overall evidence for a cross- correlation different from 
zero at the 2.7 a level, while this evidence drops to 1.5 a at z > 1.5. We focus on the capabilities 
of the ISW to constrain the behaviour of the dark energy component at high redshift both in the 
ACDM and Early Dark Energy cosmologies, when the dark energy is substantially unconstrained 
by observations. At present, the inclusion of the ISW data results in a poor improvement compared 
to the obtained constraints from other cosmological datasets. We study the capabilities of future 
high-redshift QSO survey and find that the ISW signal can improve the constraints on the most 
important cosmological parameters derived from Planck CMB data, including the high redshift dark 
energy abundance, by a factor ~ 1.5. 



I. INTRODUCTION 

The cosmic microwave background (CMB) measure- 
ments of temperature anisotropies and polarization [lj 
and the redshift-distance measurements of Type la Su- 
pcrnovae (SNIa) at z < 2 @, have established that the 
universe is undergoing an accelerated phase of expan- 
sion and that its total energy budget is dominated by a 
dark energy component. The nature of this component is 
still unknown and many observational probes have been 
proposed to test its properties and redshift evolution ei- 
ther in the standard A-Cold Dark Matter (ACDM), mod- 
ified gravity or quintessence models (for a review see Ref. 
0). While the CMB is a powerful cosmological probe 
of the universe at z ~ 1100, the anisotropies present in 
CMB data (kinetic and thermal Sunyaev-Zeldovich ef- 
fects, weak lensing and the Integrated Sachs- Wolfe effect) 
contain precious information on the large scale structure 
that formed at much lower redshift. These effects can be 
studied and detected by cross-correlating CMB data with 
tracers of the large scale structure (LSS) such as galaxies 
or quasars. Here, we will focus on the Integrated Sachs- 
Wolfe effect [|. 

The first investigations of CMB-LSS cross-correlations 
were made in Refs. @ using X-ray observations 



(HEAO data) and CMB data from the COBE satellite 
that allowed to constrain the bias of X-ray sources and 
put upper limits on the amount of a cosmological con- 
stant energy density. In a more recent series of works, a 
similar set of analyses have been carried which relied on 
CMB data from the WMAP satellite and a variety of LSS 
probes such as NVSS (NRAO VLA Sky Survey) radio 
galaxies 0,0], the Two Micron All Sky Survey (2MASS) 
H, G3, galaxies from Sloan Digital Sky Survey (SDSS) 
DR4 and DR5 [HI, [Hj]. There is overall agreement be- 
tween the different groups in finding an evidence for a 
positive ISW signal at the ~ 3 c confidence level. Even 
more recently, several groups have started to analyze the 
CMB-LSS correlations in a more comprehensive frame- 
work fully exploiting the capabilities of the SDSS (DR5, 
DR6), addressing the many possible systematic effects 
involved, modeling the statistical error bars and covari- 
ance properties of the data with different methods and 
combining the different tracers in the using Monte Carlo 
Markov Chains estimators for deriving the cosmological 
parameters (e.g. 0, 0, QJ, QJ]). 

Although the results obtained from the CMB-LSS are 
promising and have provided independent evidences for 
the presence of a dark energy component, their quan- 
titative use as cosmological probes is still to be fully 
exploited. The timely convergence of future LSS sur- 
veys like SDSS-III [11 [II] and high resolution CMB ex- 
periment like Planck 1 will offer the opportunity to fur- 
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ther explore the valuable information of the ISW and in 
general of several other LSS-CMB cross-correlations in 
a quantitative way (e.g. [HI). Present efforts are par- 
ticularly concentrated on quantifying several systematic 
effects that are possibly affecting the ISW signal such 
as noise by local variance [2fj], redshift space distortions 
[2l| . non- linear contributions [12], uncertainties in the 
bias estimates 23] and contributions due to voids and 
clusters to the overall ISW signal [24j . 

However, once the systematic effects will be under con- 
trol the ISW is likely to become a powerful cosmological 
probe able to constrain for example the contribution of 
massive neutrinos (25| or evolution of the dark energy 
component at high redshift (e.g. [11]), which is poorly 
constrained by observations at the present [I?! • In this 
paper we focus on this latter issue and push the ISW ca- 
pabilities to the highest possible redshift regime in a way 
similar to that used by Ref. jl3j, using the recently re- 
leased SDSS DR6-QSO catalog of about one million pho- 
tometrically selected QSOs in order to give constraints on 
some dark energy models. 

Among all the dynamical dark energy models, we will 
consider early dark energy (EDE) ones, in which a small 
fraction of dark energy is present up to the last scatter- 
ing surface (lss), unlike ACDM for which SIde^Iss) — 0, 
The differences between early dark energy models and 
pure ACDM are particularly evident at high redshifts, 
over a large fraction of the cosmic time, when the first 
structures form. EDE has been shown to influence the 
growth of cosmic structures (both in the linear and in 
the non-linear regime), to change the age of the uni- 
verse, to have an influence on CMB physics, to impact 
on the reionization history of the universe, to modify the 
statistics of g iant arcs in strong cluster lensing statistics 

(e.g. m, mm m m m iiisii m). 

Recently some EDE models have been investigated by 
Refs. jiH |42j, focussing in particular on using present 
and future weak lensing observables and measurements 
of the growth factors of density perturbations obtained 
via Lyman-a forest observations. 

The structure of the paper is as follows: in Sec. HTland 
Sec. IIIII we briefly review the theoretical background of 
the ISW effect and describe the quasar catalog used, re- 
spectively. Sec. IIVI and Sec. |V] contain the analysis of 
ACF of QSO number density and CCF between QSO 
number density and CMB temperature fluctuations. In 
Sec. lVII we present the theoretical framework of the early 
dark energy model and the datasets we used. Sec. IVIII 
contains the bulk of our results, while Sec. I Villi is ded- 
icated to forecasting with future datasets. We conclude 
with a discussion in Sec. IIXI 



II. ISW EFFECT 

In this section we briefly review the basics of ISW effect 
focussing on its cross-correlation with the number density 
of astrophysical sources (e.g. Refs. [HHEI). 



The temperature anisotropy due to the ISW effect is 
expressed as an integral of the time derivative of the grav- 
itational potential $ over conformal time r\ 

— (n) = -2 J h( m - V )]dT) . (I) 

For scales within the horizon, we can relate the gravi- 
tational potential $ to the comoving density field 5 by 
Poisson equation: 
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where f2 m is the ratio of the matter density to the critical 
density today, Hq is the Hubble constant today, c is the 
speed of light, z is the redshift, and k is the comoving 
wave number. From Eqs. (1 1 12)) one can appreciate that 
when a CMB photon falls into a gravitational potential 
well, it gains energy, while it loses energy when it climbs 
out of a potential well. These effects exactly cancel if the 
potential is time independent, such as the matter domi- 
nated era (S m ~ a) in which the gravitational potential 
stays constant, $ = and no ISW is produced. However, 
when dark energy or curvature become important at later 
times, the potential evolves as the photon passes through 
it. In this case, $ and additional CMB anisotropics 
will be produced. 

The ISW effect of interest here is the one produced at 
relatively late time, when the dark energy component is 
dominating the universe density budget, causing a change 
in the time dependence of the expansion rate, departing 
from pure matter dominance, and a consequent time evo- 
lution of the gravitational potentials. An early ISW, not 
considered here, is injected soon after decoupling, when 
the expansion rate time dependence is in the transition 
between radiation and matter dominance. Observing the 
late-time ISW can be a powerful way of probing dark en- 
ergy and its evolution. However, the most significant 
ISW effect contributes to the CMB anisotropies on large 
scales that are strongly affected by the cosmic variance. 
Fortunately, this problem can be solved by the cross- 
correlation between ISW temperature fluctuation and the 
density of astrophysical objects like galaxies or quasars 
(in the following calculations we will use the quasar cat- 
alog). The observed quasar (QSO) density contrast in a 
given direction fii will be: 

(5 q (rii) = / f(z)S m (&i,z)dz 



dN 

b q {z) — {z)S m (n ll z)dz, (3) 



where b q (z) is an assumed scale-independent bias factor 
relating the quasar overdensity to the mass overdensity, 
8q = b q S m , dN/dz is the normalized selection function of 
the survey. Since the density S m is related to the grav- 
itational potential the observed galaxy density will 



3 



be correlated with the ISW temperature in the nearby 
direction n 2 : 



AT 



dz 



(n 2 , z)dz 



(4) 



Given a map of CMB and QSO survey, the angular 
auto-correlation and cross-correlation functions can be 
easily expressed in the harmonic space: 
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where # = |ni — n 2 | and the auto-correlation and cross- 
correlation power spectra are given by: 



qT 



k 2 dkP(k)ll sw (k)I?(k) , (7) 



c qq - - / k l dkP(k)[i?{k)Y , 



(8) 



where P(k) is the matter power spectrum today and the 
functions // sw (fc) and J q (fc) are: 



rISW 



(k) 



d$(fc) 



ji[k\(z)}dz , 



(9) 



I?(k) = b q (z) — (z)5 m (k,z)ji[kx(z)}dz ,(10) 



where are the spherical Bessel functions, and \ is 

the comoving distance. In the following, we use the pub- 
lic package CAMB_sources 2 to calculate the theoretical 
angular auto-correlation and cross-correlation functions. 



III. QUASAR CATALOG 

We use the SDSS DR6 quasar catalog released by 
Ref. [45[ (hereafter DR6-QSO). This unique quasar cat- 
alog contains about A q i « 10 6 objects with photo- 
metric redshifts between 0.065 and 6.075, covering al- 
most all of the northern hemisphere of the galaxy plus 
three narrow stripes in the southern, for a total area of 
8417 deg 2 (~ 20% area of the whole sky). Photometri- 
cally selected QSOs have become particularly important 
in the last few years due to the higher selection efficiency 
reached, that has enabled their use for meaningful sta- 
tistical/cosmological analysis. The DR6-QSO data set 
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FIG. 1: Redshift distributions of the two quasar subsamples: 
black solid line is for the "UVX-A11" subsample while red 
dashed line is for the "All- All" subsample. We also plot the 
theoretical redshift distribution function dN/dz of two sub- 
samples: blue solid line and green dashed line are for "UVX- 
AU" and "All- All" subsample, respectively. 



extends previous similar SDSS data sets with ~ 95% ef- 
ficiency [46|, |4?J ■ The main differences are due to the fact 
that DR6-QSO probes QSOs at higher redshift and also 
contains putative QSOs flagged as to have ultra violet 
excess (UVX objects). We refer the reader to Ref. [45[ 
for a very detailed description of the object selection with 
the non-parametric Bayesian classification kernel density 
estimator (NBC-KDE) algorithm. 

For our purposes we rely on the electronically pub- 
lished table that contains only objects with the "good" 
flag with values within the range [0, 6], and we label this 
as "All-All" subsample. The higher the value, the more 
probable for the object to be a real QSO (see Section 
4.2 of Ref. [45j for details). As a further more conserva- 
tive criterion, we will present results for a subset with 
"uvxts=l", i.e. objects clearly showing a UV excess 
which should be a signature of a QSO spectrum. We 
are left with A 7 q2 w 6 x 10 5 quasars and we refer to this 
sample as "UVX-AH" subsample. This second choice is 
clearly more conservative than the first one and thereby 
results should be trusted more at a quantitative level. 

For simplicity, we assume that the redshift distribution 
dN/dz of the DR6-QSO sample is approximated by the 
function: 



dN 

dz 







■pf m+i 



(11) 



where m, (3 and z$ are three free parameters, which are 
dependent on the redshift distribution of quasar number 
density. For the "All-All" subsample, we find that m = 
2.00, f3 = 2.20, and z Q = 1.90. The mean redshift of this 
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subsample is z ~ 1.80. The "UVX-AH" subsample has 
m = 2.00, j3 = 2.20, z = 1.62 and its mean redshift 
is z ~ 1.49. The distributions are shown together with 
the redshift distribution of the quasar number density in 
Fig. Q] In this plot we do not normalize the distributions 
to be unity. 

Furthermore, we also take into account the magnifica- 
tion bias effect which could be important for the SDSS 
QSOs. In the presence of magnification bias, the relevant 
quantity entering Eq. (|10p is the function f(z) given by 

m 



dN 

/(*)=6(*)^(*) 



W{z,z')[a(z') 



1] dz^ dz ' 



(12) 

where a(z') is the slope of the number counts of the 
quasar number density as a function of flux: N{> F) oc 
F~ a . For simplicity, in our analysis we set a = 0.9 in 
the whole redshift region. Here in the flat universe the 
lensing window function W(z,z') is [lj|: 



i 



i 



X(z) x(z') 



(13) 



where x( z ) = Jo dz" / H(z") is the radial comoving dis- 
tance. 

Besides the quasar sample in the whole redshift region, 
we are also interested in the quasar subsamples in the 
low redshift z < 1.5 and the high redshift z > 1.5. In the 
"All- All" subsample, there are 7V lq i ~ 3.54 x 10 5 QSOs at 
low redshifts z < 1.5 ("All-Low") and A hq i w 6.57 x 10 5 
QSOs at z > 1.5 ("All-High"). The mean redshifts are 
0.90 and 2.28, respectively. For the "UVX-A11" subsam- 
ple, we are left with A lq2 « 2.87 x 10 5 ("UVX-Low") and 
N hq2 ~ 3.21 x 10 5 ("UVX-High") quasars. Their mean 
redshifts are 0.90 and 2.02, respectively. 

In these four cases, we also assume that the theoreti- 
cal redshift distributions are approximated by the func- 
tion Eq. (|11[) in the redshift region [21,32]. Furthermore, 
we widen these redshift distributions by two one-tailed 
Gaussians when < z < z\ and z > Z2 [47] : 



dN 

-j-{z) = zi,«2 exp 
dz 



(z - [z 1 ,z 2 ]f 
2a 2 



(14) 



where a is the dispersion between photometric and spec- 
troscopic redshifts. This technique of widening by Gaus- 
sians is consistent with determining dN / dz from spectro- 
scopic matches in the photometric redshift bin [471 ]. In 
Ref. [13] , the authors measured this dispersion in the dif- 
ferent redshift bins and found a ~ 0.2. Therefore, in our 
following calculations we simply set a = 0.2. 

We have explicitly checked that these choices for the 
parameters a and a have negligible impact on our final 
results. 




50.0 



FIG. 2: Zoom on the SDSS DR6-QSO number density map 
with iVsidc = 128 for the "All-All" subsample. This map is in 
celestial coordinates. The grid spacing is 30° and the center is 
at RA = 180 deg, Dec = 37deg, corresponding to / = 167.45 
and 6 = 75.3 degrees in Galactic coordinates. RA is increasing 
to the left. 



IV. DR6-QSO AUTO-CORRELATION 
FUNCTION 

In this section we will use DR6-QSO sample to com- 
pute the QSO auto correlation function (ACF). 



A. QSO Pixelation 

We pixelate the quasar maps using the HEALPix soft- 
ware package 3 [48j |. We use a relatively coarse resolu- 
tion: A^dc = 64, corresponding to A pix = 49, 152 pixels 
with dimensions 0.92° x 0.92°. This resolution is suffi- 
cient for the large scale correlations we are interested in. 
In FigJ^l for illustrative purposes, we show the quasar 
number density map of the "All-All" subsample in ce- 
lestial coordinates. For the lowest resolution map, there 
are only ~ 20% of these pixels actually contain quasar 
sources. 

In the low resolution map (LA S jd c = 64), there are 
many edge pixels which are only partially filled by the 
QSOs. To account for such an effect, following Ref. [HI], 
we also use a higher resolution pixelation map (HN s ide = 
512) to estimate the coverage fraction of each edge pixel. 
In order to determine the mask of the actual sky cover- 
age of the DR6 survey, we generate a random sample of 
galaxies using the DR6 database to ensure roughly uni- 
form sampling on the SDSS CasJobs website 4 . By using 
a sufficiently large number of random galaxies (between 
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13—20 million) we can make sure to have a good sampling 
when these are pixelated on the high resolution map. 

Here, we choose LN s ide = 64 and HN s i dc = 512 for the 
low and high resolution maps, respectively. We estimate 
the coverage fraction of each low resolution pixel, fi, as: 



fi = 



N. 



high 



(0 



1)4 



(15) 



where JV^fi) is the number of high resolution pixels 
within the mask for each coarse pixel i, and for choice 
made here there are 64 = (HN sidc /LN sidc ) 2 = (512/64) 2 
high resolution pixels in each coarse pixel. Using this 
method now most of edge pixels are partially covered by 
the SDSS DR6 survey, < /; < 1. 

We then correct the maps by dividing the observed 
number of quasars in each coarse pixel by the fraction 
of the sky within the pixel that was observed, yielding 
n,//,. We use the higher resolution to calculate the av- 
erage number of quasars per coarse pixel, n: 



N„ 



x64=^- 



(16) 



We also use a higher low-resolution LN s - lde = 128 and 
a higher high-resolution HN a i de — 1024 to perform all 
the calculations and find that our results are stable. 



B. ACF Estimator 

In order to measure the DR6-QSO ACF, we use the 
ACF estimator c u (9), where the index tt refers to the 
total catalog (including possible stellar contaminations): 



c"(0) 



N 



-E 



(m - fin){nj - fori) 



N, 



Em 



(17) 
(18) 



where fi is the pixel coverage fraction, rii is the number 
of quasar sources in each pixel, and fi is the expectation 
value for the number of objects in the pixel. The sum 
runs over all the pixels with a given angular separation. 
For each angular bin centered around 9, Ng is the num- 
ber of pixels pairs separated by an angle within the bin, 
weighted with the coverage fractions. 

This estimator is equivalent to the one used in Ref . [49[ 
in which the authors construct the ACF from the counts 
of data-data, random-data and random-random number 
density pairs: 



c u (9) 



QQ{9) + RR{9) - 2QR{9) 
RR{9) 



(19) 



where Q and R denote the data point and random point, 
respectively. In the following calculations, we also use 
this estimator for cross-checking and find consistent re- 
sults. 



Because we pixelate the quasar map using a low reso- 
lution iVgido = 64, in which the pixel size is 55' [Hj], we 
use Ah = 12 angular bins in the range 1° < 9 < 12° and 
a linear binning in our calculation. The choice of binning 
does not affect the results significantly. 



C. Covariance Estimator 

We estimate the covariance matrix of the data points 
using jackknife resampling method [5fj| . This method is 
to divide the data into M patches, then create M sub- 
samples by neglecting each patch in turn. These patches 
have roughly equal area. In practice, we firstly list the 
whole set of pixels covered by the survey, and then divide 
them into M = 30 patches that do not have very similar 
shape, but have roughly equal area (i.e. equal number of 
pixels) . 

The covariance estimator reads: 



M - 1 ^ 



M 



k=l 



C* 



C t k t (9 ] )-C tt {9 3 



where 



(20) 

) are the observed ACF of the M subsamples 
in the i-th angular bin and C tt (9i) are the mean ACF over 
M realizations. The diagonal part of these matrices gives 
the variance of the ACF in each bin = of, while the 
off-diagonal part represents the covariance between the 
angular bins. We also change the number of patches M 
and verify that the covariance matrix is stable. We refer 
the readers to Ref. [lj| for a more extensive comparison 
between several covariance estimators. 



D. Stellar Contamination 

Although the SDSS DR6-QSO catalog we use has very 
high efficiency in the selection algorithm, stars are point- 
like sources that inevitably contaminate the catalog. An 
estimate of the level of stellar contamination can be com- 
puted using the fact that the correlations properties of 
stars are very different from those of QSOs. In particu- 
lar, a nearly flat contribution, up to large angular scales 
(that correspond to small physical Galactic distances) is 
expected from stars of our own galaxy. 

If we substitute rn — > anj + (1 — <z)nf, where a is the 
efficiency of the quasar catalog, the ACF estimator of 
Eq. (fT7|) becomes: 



z u {9) 



a 2 c"(0) + (1 - afc ss {9) + e{6) 



(21) 



where c qq {9) and c ss (9) are the intrinsic ACF of QSOs 
and stars, respectively, and e(9) is a tiny offset arising 
from cross-terms (47l |: 

(m 0{ 2 JQS(9) + RR(9)-QR(9)-SR(9)\ 
e(9) = 2 (a — a ) 1 ' 



2(a-a 2 ) 



RR{9) 
f l n)(n s ] - fjfi) 
fi 2 



(22) 
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FIG. 3: The auto-correlation function of star catalog c ss (8). 
And the errors are also derived from the jackknife method. 



FIG. 4: The auto-correlation function of quasar "UVX-A11" 
subsample measured for all the sample, for a single foreground 
mask and for all masks joint. No error bars are reported. 



If the efficiency of the catalog is high enough, i.e. a > 
90%, this cross-term e{9) should be close to zero and can 
be neglected safely in the analysis. Otherwise, e(9) will 
become large and is comparable to the intrinsic ACF of 
stars c ss (9). Thus, in order to correctly determine the 
efficiency of the catalog with a large stellar contamina- 
tion, such as the "All- All" subsample, we have to take 
this term into account. We will show this effect in the 
following sections. 

To compute the stellar contamination, we extract a 
large number (~ 8 x 10 4 ) of stars in SDSS DR6 sur- 
vey in the magnitude range 16.9 < g < 17.1 from the 
CasJobs website and compute the ACF of stars c ss (9) us- 
ing Eq. (TT7) . In Fig.[3]we show the stellar ACF and its er- 
rors which are also measured with the jackknife method. 
We can find that the stellar contamination will domi- 
nate the ACF of the total catalog at the largest scales 
where the ACF of QSOs c qq (6>) -> 0. The star contri- 
bution is indeed not perfectly flat and retains a small 
slope even at large angles, dropping from c ss (l°) ~ 0.3 
to c ss (T2°) ~ 0.15, which is consistent with the result 
of Ref. [43| , showing that probably our star catalog is 
also contaminated at some level by QSOs. However, the 
contribution at large angular scales could be robustly es- 
timated and removed. 



E. Systematic Errors 

Several systematics effects, including the galactic ex- 
tinction by dust, sky brightness, number of point sources 
and poor seeing, could potentially affect both the ob- 
served ACF and CCF. These systematics are investigated 
in detail in Refs. [HI, [TH, H3|. We also checked for their 
contribution in our calculation and in particular we con- 
sider extinction and point sources contamination, that 



are believed to affect most the measured correlations. 

Basically we recompute the ACF by considering only 
objects within the pixels with A g < 0.18 (the g-band 
Galactic extinction) for the reddening, while for the point 
sources we remove the pixels with more than twice av- 
erage number density of sources, rii < 2 h. These masks 
will remove about ~ 20% of the considered area. Among 
these two systematics we find that the extinction have 
the major effect on the derived ACF and we plot the 
observed ACF with and without these foreground masks 
for the "UVX-Air subsample in FigOU We also plot the 
effect of extinction on the ACF for the "UVX-Low" and 
"UVX-High" subsamplcs in Fig[5] Therefore, we will use 
the extinction mask A s < 0.18 only in our following anal- 
ysis. Our results are in quantitative agreement with those 
of Ref. [IH based on SDSS DR4 quasar catalog. 



F. Auto-Correlation Functions 

In this subsection, we summarize the whole calcula- 
tions of DR6-QSO ACF. For our purpose we construct 
six different quasar subsamples from the SDSS DR6- 
QSO catalog. Firstly we pixelate the quasar maps using 
HEALPix software and remove the ~ 20% highest con- 
taminated pixels with A g > 0.18. Next we use the ACF 
estimator Eq. (jTTJ) and the covariance estimator Eq. (j2"D|) 
to calculate the ACF and its covariance matrix for each 
quasar subsample. Finally we consider the contribution 
of stellar contamination and also the cross terms when 
the efficiency is not high enough. Note that in the fol- 
lowing analysis we take in account the window function 
w(9) associated with our pixelation. 

In Fig. [5] we plot the observed ACF for these six differ- 
ent quasar subsamples. We find that for the "UVX-All" 
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FIG. 5: The auto-correlation function of quasar "UVX-Low" 
and "UVX-High" subsamples measured for the extinction 
mask. No error bars are reported. 



subsample the values of ACF are consistent with previous 
works [13l. Il4j. However, if we consider the high redshift 
subsample, the observed ACF becomes large. The rea- 
son is that in this case the subsample has a higher mean 
redshift. At this high mean redshift, the bias and stellar 
contamination will become larger and produce a larger 
ACF (the efficiency is smaller). Moreover, due to their 
higher mean redshift, larger bias and smaller efficiency, 
three "All-x x x" subsamples also have larger observed 
ACFs than the corresponding three "UVX-x x x" sub- 
samples. In Figj6]we show the theoretical ACFs for differ- 
ent cases. Obviously all the plots depend on the chosen 
cosmology. Here we use the best fit model of WMAP 
five-year data Q: n h h 2 = 0.02267, n c h 2 = 0.1131, 
r = 0.084, h = 0.705, n s = 0.96 and A s = 2.15 x 10~ 9 at 
k = 0.05 Mpc -1 . For the other two free parameters, bias 
and efficiency, we choose the best fit values of Table II, 
obtained from our calculations in the ACDM framework 



that will be shown in the following sections. We find 
that there is good agreement between the theory from 
the WMAP5 best fit model and the observed ACF in 
each quasar subsample. 



QSO-CMB CROSS-CORRELATION 
FUNCTION 



For the cross-correlation analysis, we use the WMAP 
Internal Linear Combination (ILC) map derived from the 
five-year WMAP data, with HN B id c = 512 provided by 
the WMAP team [5l|, shown in Fig|71 in the same re- 
gion as the SDSS DR6-QSO data. This ILC map was 
already built to minimize the Galactic and other fore- 
ground contaminations. In our calculations, we down- 
grade the ILC map to the low resolution LN^p, = 64. 
For the WMAP mask, we use the "KQ75" mask [52] cor- 
responding roughly to the "KpO" cut in the three-year 
data release. This mask map has also N s ia c = 512. We 
down-grade this map to the low resolution and set the 
weight w T — for all pixels including at least one masked 
high resolution pixel [7j ■ 

To measure the CCF between the SDSS DR6-QSO 
number density map and the WMAP ILC map, we rely 
on the following estimator: 



(0) 



V(Ti-T) 3 _ 3 

— * n 



(23) 



1,3 



where Ti is the CMB temperature in the i-th pixel and T 
is the mean (monopole) value for the CMB temperature 
in the unmasked area. N$ is the number of pixels pairs, 
which has been defined in the previous Section. We also 
use Nb = 12 angular bins in the range 1° < 8 < 12° 
and a linear binning in our calculation. We also check 
the contributions from the extinction and point sources 
contamination in Fig[9l and find that the extinction have 
the major effect on the CCF. We also plot the effect of 
extinction on the CCF for the "UVX-Low" and "UVX- 
High" subsamples in FigdU] Therefore, we apply the 
reddening mask with A g < 0.18 to be consistent with the 
ACF measurements. 

We also consider the possible contributions from stars 
in the CCF calculations. The stellar contamination has 
to be subtracted from the total CCF and the QSO- 
temperature correlation that will be compared with the 
theoretical models becomes: 



(0) 



c tx (0)-(l-a)c si (0) 



(24) 



where the star-temperature cross-correlation has to be 
also estimated. 

The covariance matrices are also calculated by using 
jackknife resampling method, Eq. (f2"D|) . We firstly list all 
the pixels covered by the survey in the quasar map, and 
divide them into M — 30 patches. Then we create M 
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FIG. 6: The observed auto-correlation function (ACF) of different quasar subsamples. Top three panels are for "UVX-A11", 
"UVX-Low" and "UVX-High" subsamples. Bottom three panels are for "All- All", "All-Low" and "All-High" subsamples. All 
errors are estimated with jackknife method. The red solid lines are the theoretical predictions from the best fit model of 
WMAP five-year data, while the black dashed lines and the blue dash-dot lines are its component (theoretical quasars ACF 
and stellar contamination). We also plot the contributions from the cross term e(8) (cyan short dashed lines) in the "All- All" 
and "All-High" subsamples in which the efficiencies become small. 




FIG. 7: WMAP Internal Linear Combination map with 
■Aside — 128 in celestial coordinates zoomed on the same re- 
gion probed by SDSS DR6-QSOs. The grid spacing and center 
are the same as in Fig(2] 



subsamples by neglecting each patch in turn and discard 
the same pixels in the WMAP ILC map. 

In Fig. [5] we plot the observed CCF for these six dif- 
ferent quasar subsamples, together with the theoretical 



CCF predictions of WMAP5 best fit cosmological model. 
There is good agreement between the theoretical and ob- 
served CCFs. 



The values of CCF of the "UVX-AH" subsample are 
consistent with previous works [HI, llij ]. We find 0.22 ± 
0.10/jK and 0.26 ± 0.09^K for the "UVX-A11" and "All- 
All" subsamples at 1 deg, respectively. The median value 
of "All- All" subsample is larger than one of "UVX-AH" , 
due to the larger bias and stellar contamination. The 
larger number of quasars in "All- All" subsample allow 
to slightly shrink the error bars. In particular, we find 
that the high redshift subsample z > 1.5 also gives a 
non-vanishing cross-correlation signal, which is consis- 
tent with the results of Ref. [l5j]. One of the possible 
reasons could be that the high redshift subsample has 
a larger bias and stellar contamination. Another specu- 
lative possibility is also that a larger value for the ISW 
CCF than that implied by ACDM could reflect a different 
underlying cosmological model, characterized by an early 
departure from matter dominance and onset of acceler- 
ation, as it happens in early dark energy and modified 
gravity scenarios. 
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FIG. 8: The observed cross-correlation function (CCF) between different quasar subsamples and WMAP ILC map. Top three 
panels are for "UVX-AH", "UVX-Low" and "UVX-High" quasar subsamples. Bottom three panels are for "All- All", "All-Low" 
and "All-High" quasar subsamples. All errors are jackknife estimated. The red solid lines are the theoretical predictions from 
the best fit model of WMAP five-year data. 



VI. METHOD AND DATASETS 



o 

u 0.1 



■ 


1 1 








■ 


1 1 1 1 

No Mask 


• 










• 


Point Sources 












A 


Extinction 




1 








T 


All Mask 




■ 




























♦ 


1 














* 


1 


s 


• 










+ 




♦ . 
















: • 














♦ * • . 



6 (deg) 



FIG. 9: The cross-correlation function of quasar "UVX-AU" 
subsample measured for all the sample, for a single foreground 
mask and for all masks joint. No error bars are reported. 



A. Parametrization of Early Dark Energy 

Due to the possible non-zero ISW signal at high red- 
shift provided by the early dark energy model, we will 
also consider here the EDE mocker model introduced in 
Ref. [H: 



wede(o) 



1 + w 



(25) 



where a = l/(l + z)is the scale factor, wq is the present 
equation-of-state of dark energy and C characterizes the 
"running" of the equation of state. Consequently, the 
evolution of dark energy density can easily be obtained 
via energy conservation as: 



Pede(q) 
Pede(1) 



= [(1 



w )a 



-c 



w 



,3/C 



(26) 



In Fig. [IT] we plot the dark energy density (upper part 
of the panel) and equation of state (bottom part of the 
panel) as a function of rcdshift for two different models: 
pure ACDM (blue lines) and EDE (red lines) which has 
(w ,C) = (-0.95,2.50). All these models fit the CMB 
and the lower redshift SNIa constraints very well [4l| . 
In the panel we also show as a shaded vertical band the 
region in the redshift range z = — 5. This is the redshift 
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FIG. 10: The cross-correlation function of quasar "UVX- 
Low" and "UVX-High" subsamples measured for the extinc- 
tion mask. No error bars are reported. 



FIG. 11: The evolution of dark energy density and equation 
of state for two models, wo = —0.95 and C = 2.5 (EDE) and 
pure ACDM model. The shaded area represents the redshifts 
that is probed by SDSS DR6-QSOs (z ~ - 5). 
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range we will be (mainly) focussing on in the rest of the 
paper. 

We stress that this is just one of the possible parame- 
terizations for early dark energy models, another one has 
been suggested by Ref. [53| and recently used in Ref. • 
However, we prefer to use the mocker model in order to 
compare with the findings of Ref. and because this 
parametrization has a smooth redshift derivative at low z 
for w(z). We note that one of the most important param- 
eters is the amount of dark energy during the structure 
formation period and this is given by (a oq is the matter- 
radiation equality scale factor): 

,o 

^EDE, s f = -(In a cq ) / Oede (a) din a , (27) 

J In a cq 

and we will also quote this value in the rest of the paper 
in order to com par e with other works and constraints as 
well (e.g. Refs. (HIH). 



FIG. 12: The late-time ISW temperature power spectrum in 
two different dark energy models: wo — —0.95 and C = 2.5 
(EDE) and pure ACDM model. 



In Fig. [T2] we plot the late-time ISW temperature 
power spectrum in the two different dark energy models 
analyzed. We find that the early dark energy model pro- 
duces a larger contribution to the CMB primary power 
spectrum at low multipoles than the ACDM model due 
to the different low redshift evolution: differences up to 
a factor ~ 2 are present for I values smaller than 10, 
reflecting the fact that ISW cross-correlation is in place 
earlier than in ACDM, due to the early departure from 
matter dominated expansion rate. 

We also focus on the ISW effect at different redshifts 
investigating EDE models. We modified the CAMB code 
[HI to calculate the contributions of ISW effect to the 
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FIG. 13: Top: Ratio of the CMB temperature power spectrum with and without the ISW effect at different low redshifts in two 
different dark energy models: EDE (red dashed lines) and ACDM (black solid lines). Bottom: The cross-correlation functions 
C aT {9) between EDE Mocker model (red dashed lines) and pure ACDM model (black solid lines). The redshift distribution 
dN/dz is assumed to be a Gaussian function centered at three different mean redshifts zo = 1, 2, 3 comparison with a a = 0.5. 
We fix bias to b — 2.2 for illustrative purpose. 



power spectrum prior to some redshifts z*, and neglect 
the ISW effect at low redshifts: 



AT 



isw 



(n) = , when z < z* 



(28) 



In Fig. [T3] we plot the ratio of the CMB temperature 
power spectra with and without the ISW effect at differ- 
ent low redshifts and in two different dark energy models. 
We find that at the lowest redshift z < 1, the ISW contri- 
bution of ACDM is larger than the EDE model. However, 
when the redshift becomes larger, the EDE model con- 
tributes more to the ISW effect and to the CMB temper- 
ature power spectrum. Therefore, when calculating the 
cross-correlation power spectra between ISW and QSO 
survey with different mean redshifts, we can see that the 
higher the mean redshift of QSO survey is, the larger the 
cross-correlation C ; qT becomes, which is shown in Fig. 1131 
From these plots it is clear how the cross-correlation sig- 
nal between QSO number density and CMB temperature 
differs in the two models up to a factor three (the dif- 
ferences at higher redshifts z > 3 becoming even larger 
than those shown here). Moreover, the redshift evolution 
of the ISW signal is quite different in the two models and 
this is promising for future studies that aim at investi- 
gating high-redshift deviations from the standard model. 



B. ISW Likelihood Function 

In Refs. [Hill, the authors firstly used the QSO ACF 
data only to determine the bias b, when fixing the other 
parameters and then they used the CMB-QSO CCF data 
to constrain the significance of ISW signal or other cos- 
mological parameters. However, in order to make the 
whole analysis consistent, in our calculations we use both 
the ACF and CCF data to constrain all the parameters, 
including the bias b, efficiency a, ISW amplitude A amp , 
which is defined in Eq. (|3"3"j). and other cosmological pa- 
rameters, simultaneously. 

We then compare the theoretical ACF c qq (8) and 
CCF c qT (0) with the observed values of ACF c qq (6) and 
CCF c qT (0), respectively, through the Gaussian likeli- 
hood function: 

£acf = (2n)- N / 2 [detC t3 ]-^ 2 



x exp 



v ^ i (gr-cr)(s qq -s qq ) 



N/2, 



£ccf = (27r)-^[detCy- 1/2 



x exp 



2-^ 9 



,(29) 



(30) 



hJ 



where Cy and C'^ are the DR6-QSO auto-correlation 
function and CMB-QSO cross-correlation function co- 
variance matrix. 
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C. Other Current Datasets 

Besides the ACF and CCF data, we will rely here on 
the following cosmological probes: i) CMB anisotropics 
and polarization; ii) baryonic acoustic oscillations in the 
galaxy power spectra; iii) SNIa distance moduli. 

In the computation of CMB power spectra we have 
included the WMAP five-year (WMAP5) temperature 
and polarization power spectra with the routines for 
computin g th e likelihood supplied by the WMAP team 

BAOs (Baryonic Acoustic Oscillations) have been de- 
tected in the current galaxy redshift survey data from the 
SDSS and the Two-de gree Field Galaxy Redshift Survey 
(2dFGRS) [H, S3, W&m- The BAO can directly mea- 
sure not only the angular diameter distance, Da{z), but 
also the expansion rate of the universe, H(z), which is 
powerful for studying dark energy [63| . Since current 
BAO data are not accurate enough for extracting the in- 
formation of Da{z) and H(z) separately 
only determine an effective distance [59| : 



one can 



D v (z) 



{X + zfD\{z) 



H{z) 



1/3 



(31) 



In this paper we use the Gaussian priors on the distance 
ratios r s (z d ) / D v (z): 

r s {z d )/D v (z = 0.20) = 0.1980 ± 0.0058 , 
r s (z d )/D v {z = QM) = 0.1094 ± 0.0033 , (32) 

with a correlation coefficient of 0.39, extracted from the 
SDSS and 2dFGRS surveys [12], where r s is the comoving 
sound horizon size and Zd is drag epoch at which baryons 
were released from photons given by Ref. [65| . 

The SNIa data provide the luminosity distance as a 
function of redshift which is also a very powerful mea- 
surement of dark energy evolution. The supernovae data 
we use in this paper are the recently released Union com- 
pilation (307 samples) from the Supernova Cosmology 
project [2], which include the recent samples of SNIa from 
the (Supernovae Legacy Survey) SNLS and ESSENCE 
survey, as well as some older data sets, and span the 
redshift range < z % 1.55. In the calculation of the 
likelihood from SNIa we have marginalized over the nui- 
sance parameter as done in Ref. 66]. 

Furthermore, we make use of the Hubble Space Tele- 
scope (HST) measurement of the Hubble parameter 
Hq = 100 h km s _1 Mpc -1 by a Gaussian likelihood 
function centered around h = 0.72 and with a standard 
deviation a = 0.08 [671 ]. 



D. Future Datasets 



TABLE I. Assumed experimental specifications for the mock 
Planck-like measurements. The noise parameters At and 
Ap are given in units of /^K-arcmin. 



/sky 


Z max (GHz) 


^fwhm 


A T 


Ap 


0.65 


2500 100 


9.5' 


6.8 


10.9 




143 


7.1' 


6.0 


11.4 




217 


5.0' 


13.1 


26.7 



For the simulation with Planck [68j . we follow the 
method given in Ref. [69| and mock the CMB temper- 
ature (TT) and polarization (EE) power spectra and 
temperature-polarization cross-correlation (TE) by as- 
suming a given fiducial cosmological model. In Table 
I, we list the assumed experimental specifications of the 
future (mock) Planck measurement. 

The proposed satellite SNAP 5 (Supernova / Acceler- 
ation Probe) will be a space based telescope with a one 
square degree field of view that will survey the whole 
sky. It aims at increasing the discovery rate of SNIa to 
about 2000 per year in the redshift range 0.2 < z < 1.7. 
In this paper we simulate about 2000 SNIa according 
to the forecast distribution of the SNAP [70[ . For the 
error, we follow the Ref. [70j which takes the mag- 
nitude dispersion to be 0.15 and the systematic error 
CT sys = 0.02 x z/1.7. The whole error for each data is 

given by CT mag (z l ) = ija^jzi) + 0.15 2 /m , where m is 
the number of supernovae of the z'th redshift bin. 

For the future ISW ACF and CCF data, we simulate 
two mock datasets from the best fit values of data combi- 
nation WMAP5+BAO+SNIa+"UVX-All" in the ACDM 
and EDE models, respectively. We also set b = 2.0, 
a = 97% and A amp = 1 with error bars on these param- 
eters that are reduced by a factor three. We use directly 
the covariance matrix taken from present data and divide 
it by a factor nine. This improvement could be achiev- 
able by next generation of large-scale surveys such as 
the Large Synoptic Survey Telescope (LSST, 7l]), the 
Panoramic Survey Telescope and Rapid Response Sys- 
tem (Pan-STARRS, [72j]) and the Dark Energy Survey 
(DES; The Dark Energy Survey Collaboration 2005) are 
likely to allow for an order of magnitude improvement 
in the number of QSOs which will account for the factor 
three considered here. Also, the coverage of the sky frac- 
tion is at present of the order of 20% by SDSS DR6, so 
an all-sky survey will already give a factor two improve- 
ment (we assume that the error bars scale as ^/Aqsos, 
which is reasonable if these are independent). Another 
possible improvement in this direction is the use of type-2 
QSOs instead of the type-1 used here that could further 
decrease the error bars of the sample (see the discussion 
in H|). 



In order to forecast future measurements we will use 
the same observables as before without BAO. 



5 Available at http://snap.lbl.gov/ 
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FIG. 14: Two-dimensional (marginalized) contours of flmo 
and constant Wdb from the ISW data only. 



VII. NUMERICAL RESULTS 



In our analysis, we perform a global fitting using 
the CosmoMC package [73[ a Monte Carlo Markov Chain 
(MCMC) code, which has been modified to calculate the 
theoretical ACF and CCF. We assume purely adiabatic 
initial conditions and a flat universe, with no tensor con- 
tribution. We vary the following cosmological parame- 
ters with top-hat priors: the dark matter energy density 
fl c h 2 G [0.01,0.99], the baryon energy density f^bft 2 G 
[0.005,0.1], the primordial spectral index n s G [0.5, 1.5], 
the primordial amplitude log[10 10 ^ s ] G [2.7,4.0] and the 
angular diameter of the sound horizon at last scattering 
6 G [0.5, 10]. For the pivot scale we set k sQ = 0.05 Mpc -1 . 
When CMB data are included, we also vary the optical 
depth to reionization r G [0.01,0.8]. We do not consider 
any massive neutrino contribution. From the parame- 
ters above the MCMC code derives the reduced Hub- 
ble parameter Hq, the present matter fraction il m o, erg, 
and f2EDE,sf, so i these parameters have non-flat priors 
and the corresponding bounds must be interpreted with 
some care. There are three more parameters related to 
the ACF and CCF data: the constant bias b, the effi- 
ciency of quasar catalog a and the ISW amplitude A amp 
which is defined by: 



<? qT (0) = 4 unp C qT (0) 



(33) 



where C lqT and C qT are the observed and theoretical 
CCF. In addition, CosmoMC imposes a weak prior on the 
Hubble parameter: h G [0.4,1.0]. Furthermore, we also 
include the perturbations of dynamical dark energy mod- 
els generally as done in Refs. [74l. l75l. |76| . 



TABLE II. The results of quasar bias and efficiency of 
quasar catalog in two dark energy models. 



Subsample 


Mean Redshift 


z Bias b 


Efficiency a 




ACDM model 




All-All 


1.80 


3.79 ± 0.34 


88.1% ±0.5% 


All-Low 


0.90 


1.46 ±0.29 


97.1% ±0.5% 


All-High 


2.28 


5.01 ±0.81 


80.9% ± 0.7% 


UVX-A11 


1.49 


2.18 ±0.22 


96.8% ± 0.5% 


UVX-Low 


0.90 


0.92 ±0.30 


97.5% ± 0.8% 


UVX-High 


2.02 


2.87 ±0.39 


95.1% ±0.6% 




EDE Mocker model 




UVX-A11 


1.49 


2.33 ±0.31 


96.8% ± 0.5% 


UVX-Low 


0.90 


1.05 ±0.33 


97.5% ± 0.8% 


UVX-High 


2.02 


2.98 ± 0.43 


95.1% ±0.6% 



A. ISW Only 

Firstly, we use the ACF and CCF data of "UVX-A11" 
subsample only to constrain the dark energy model with 
a constant equation of state u>de- Because at present 
the ISW ACF and CCF data have relatively large error 
bars that do not allow to use this measurement in a very 
competitive way compared to other cosmological probes, 
in our calculation we have fixed the other cosmological 
parameters to the WMAP5 best fit values, as well as the 
bias, efficiency and amplitude parameters. In Fig. [TJ] 
we plot the two-dimensional constraints on (wDE,^mo)- 
The results show that the ISW data only could give very 
weak constraints on the background parameters: u>de = 
— 1.24 ± 0.62 (la) and the 2 a upper limit of il m o < 
0.55, which is consistent with the pure ACDM model. 
This result is also consistent with some previous works 



B. Quasar Bias and Efficiency 

In our calculations we simply assume that the quasar 
bias is constant in the redshift region considered. Table 
II shows the results of quasar bias and efficiency of six 
subsamples in two dark energy models. 

For the "UVX-A11" subsample, we obtain a value for 
the bias b = 2.18 ± 0.22 (1 a) in the ACDM framework, 
also consistent with previous measurements [H, H3, [zl, 
[80j . However, we find that this bias result depends on 
the dark energy model we choose. In the EDE Mocker 
model, the bias becomes b = 2.33 ± 0.31 (lcr). Thus, 
the mean value of quasar bias becomes slightly higher in 
an EDE cosmology by ~ 10%. This effect is generally 
present in the subsamples at different redshifts. 

On the other hand, we find that the efficiency of "UVX- 
All" subsample is a — 96.8 ± 0.5%, also consistent with 
Ref. [45| . This efficiency does not depend on the different 
dark energy model, which is reasonable since this is an 
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TABLE III. Estimates of the quasar bias, from five 
photometric redshift bins. 



Redshift Bins 


Mean Redshift z 


Bias b 


0.75 <z< 1.25 


1.02 


1.58 ±0.39 


1.25 < z < 1.55 


1.43 


1.85 ±0.78 


1.55 < z < 1.95 


1.73 


2.23 ±0.87 


1.95 < z < 2.20 


2.06 


3.28 ± 0.86 


2.20 < z < 4.00 


2.46 


5.29 ± 1.31 




intrinsic feature of the quasar catalog. 

In the high redshift subsample, the quasar bias will also 
become higher and the efficiency becomes clearly smaller. 
As we mentioned before, at high redshift the stellar con- 
tamination will become large. The cross term e(9) will be 
comparable to the star ACF c ss (9) in Eq. (j21"j) . In order 
to obtain the correct efficiency, we should consider this 
cross terms in our calculations. If we neglect this cross 
term, the efficiency will be incorrectly suppressed, since 
in this case one should need a larger stellar contamination 
to contribute more ACF. For example, when using "All- 
All" subsample to do the calculations, if we neglect the 
cross terms e(fl), the efficiency is 82.5%, which is consis- 
tent with Ref. [43]. However, if taking e(9) into account, 
the efficiency rises to 88.1%. 

From the above analysis, we find that the quasar bias 
could evolve with redshift. Models in which the bias 
varies with redshift have been proposed and investigated 
(e.g. Refs. [8lL l82l]). Here, we also give the constraint 
on the evolution of the QSO bias with redshift. In our 
calculation we are not using more sophisticated models 
as those in Refs. [H, [H, [8J]. We only split the whole 
quasar "UVX-A11" subsample into five photometric red- 
shift bins: [0.75,1.25], [1.25,1.55], [1.55,1.95], [1.95,2,2] 
and [2.2,4.0], containing ~ 10 5 objects each, and assume 
a constant bias in each of them. Then we calculate their 
auto correlation function to determine their quasar bias. 
In this analysis we do not include data at z < 0.75, since 
in this low redshift the subsample is highly contaminated 
by the galaxies [47l[85|. 

In Table III and Fig. [TS] we show the estimates of the 
quasar bias in these five photometric redshift bins. We 
find that our results are in overall agreement with previ- 
ous works (47l [8fj|. Wi\ . In Fig. [15] we also plot the empir- 
ical fit derived by directly measuring the real-space clus- 
tering of spectroscopically confirmed QSOs in Ref. [7{J. 
Our results are consistent with this empirical fit very 
well. Furthermore, this analysis also allow to estimate 
important quantities such as the mass of the host haloes, 
QSO duty cycle and the Mass-Luminosity relation (see 
Ref. H). 
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FIG. 15: Estimates of the quasar bias and their 1 a error 
bars from five photometric redshift bins. The red solid line 
b(z) = 0.53 + 0.289(1 ± zf is the empirical fit derived by 
Ref. 0. 



C. Amplitude of ISW Signal 

In previous works [TJ, Il5| . cosmological parameters 
were fixed and the amplitude of ISW signal was cal- 
culated from the cross-correlation data only. However, 
generally speaking this amplitude parameter is highly 
dependent on the cosmological parameters we use. Cur- 
rently, the constraints on the cosmological parameters 
are not affected significantly by fixing the three param- 
eters; however, this is because the present data are not 
very constraining. On the other hand, in our calculations 
we find that the constraints on bias and amplitude are 
different between LCDM and EDE models (see Table II 
and IV). This fact demonstrates that keeping these pa- 
rameters fixed will lead to biased results. In other words, 
had we used the bias and amplitude obtained from the 
LCDM model to constrain the EDE model, we would 
have obtained different, biased results. Therefore, here 
we constrain all these parameters at the same time using 
the MCMC. 

The results are shown in Table IV. The significance 
in sigmas is obtained by the simple calculation, S/N = 
A/ a a- For the "UVX-AH" subsample the significance of 
ISW is about 2.7 a, which is consistent with other works 
[l3l Ibij . Furthermore, we also find that the low redshift 
subsample also gives about 2.3 a significance, which is 
reduced as compared to the "All" sample due the smaller 
quasar number of the low redshift subsample. 

The most interesting result is that the high redshift 
z > 1.5 subsample reports a signal of ISW effect at 
~ 1.5 a confidence level. As we know, in the pure ACDM 
model with the bias b ~ 1, the prediction of ISW effect 
should be close to zero in the high redshift ^4 amp ~ 0. 
However, in our analysis we obtain a higher mean value 
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TABLE IV. The results of amplitude of ISW signal and 
their significance in two dark energy models. 



Subsample 


Amplitude A am p 


S/N 




ACDM Model 




All-All 


1.51 ±0.73 


2.1(7 


All-Low 


3.05 ± 1.31 


2.3 a 


All-High 


2.32 ± 1.53 


1.5(7 


UVX-A11 


2.06 ± 0.75 


2.7a 


UVX-Low 


4.36 ± 1.82 


2.4 a 


UVX-High 


2.49 ± 1.57 


1.6a 


EDE Mocker Model 




UVX-A11 


1.91 ±0.75 


2.5 a 


UVX-Low 


3.96 ± 1.78 


2.2 a 


UVX-High 


2.22 ± 1.47 


1.5a 




0.2 0.25 0.3 0.35 65 70 75 80 




0.2 0.25 0.3 0.35 0.7 0.75 0.8 0.85 0.9 



and a smaller error bar. Consequently, the amplitude 
of ISW signal is larger than a null-detection at about 
~ 1.5<7 confidence level. As we discussed before, this in- 
teresting non-vanishing ISW signal in the high redshift 
subsample may be due to the larger bias and stellar con- 
tamination at this high redshift, or some different un- 
derlying cosmological models reflecting an early depar- 
ture for matter dominance in the expansion rate. The 
large number of quasar in the sample is also helpful in 
shrinking the error bar and enhancing the significance of 
the signal. Another intriguing result that we can appre- 
ciate from Table IV is that the discrepancies with the 
value A amp = 1, which labels the perfect agreement be- 
tween theory and data, are more severe for the lower 
redshift ranges than for those at high redshift. All the 
high-redshift samples are in fact in agreement at the 1 a 
level with the theoretical predictions either for the EDE 
or ACDM cosmologies. As for the low-redshift samples 
the discrepancies are at the ~ 2 a level confirming the 
results of Ref. [Hj]. 

As we show in Fig. [131 at high redshift the QSO- 
temperature cross-correlations of the EDE model will 
be larger than that predicted in the ACDM framework. 
When combining the obtained value for the bias, the 
EDE model should give larger CCF values than those 
of ACDM model for a given cosmology. Thus, in our cal- 
culations we find that the amplitude of the ISW signal 
in EDE model, inversely proportional to the theoretical 
CCF values (see Eq. (|33|) L is slightly smaller than the 
one of ACDM model. Although other mechanisms could 
be in place to explain this discrepancy, either involving 
not properly understood systematic effects or large scale 
structures such as super-clusters or large voids [24| , a not 
negligible amount of dark energy at high redshift could 
also help in reducing the statistical significance of this 
result. 



FIG. 16: Marginalized one-dimensional and two-dimensional 
likelihood (1,2a contours) constraints on the parameters 
fi m , Ho and ag in the ACDM model from different 
present data combinations: WMAP5 only (black solid lines), 
WMAP5±ISW (red dashed lines) and All datasets (blue dash- 
dot lines). 



D. Cosmological Constraints 

Finally we present constraints on the cosmological 
parameters from the ISW data, combining with the 
WMAP5, BAO and SNIa datasets, in two dark energy 
models. Here we only use the ACF and CCF data from 
the "UVX-A11" subsample. The other five subsamplcs 
also give similar results. In our calculations we do not 
follow the previous works [TUHH by fixing the other three 
parameters, b, a and j4 a mp, to be their best fit values. In 
Table V we show the constraints on some related cos- 
mological parameters from three different data combina- 
tions: WMAP5 only, WMAP5±ISW, and All datasets. 
And we particularly pay attention to the ISW contri- 
bution by comparing the results between WMAP5 and 
WMAP5±ISW. 



1. AC 'DM Model 

Firstly, we consider the pure ACDM model. In Fig. [TBI 
we show the one dimensional marginalized likelihood dis- 
tributions of some cosmological parameters from three 
data combinations. From Table V we can find that 
the combined constraints from WMAP5±ISW are only 
slightly improved over using WMAP5 only, since at 
present constraints from the ISW data are still very weak 
and in the calculations we only consider the quasar cata- 
log and neglect other low redshift tracers which could give 
powerful ISW constraints 0, EH EH • We also show the 
two dimensional contour in the (fi m ,(7g) panel. When 
adding the ISW data, the constraint becomes slightly 
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TABLE V. Constraints on the ACDM and early dark energy model from the current observations. Here we show the mean 
values and la error bars. For some parameters that are only weakly constrained we quote the 95% upper limit. 



Parameter WMAP5 Only WMAP5+ISW All Datasets 

ACDM Model 

n m 0.261 ±0.030 0.261 ±0.028 0.273 ± 0.019 

a s 0.797 ±0.035 0.795 ± 0.032 0.805 ± 0.027 

Ho 71.4 ±2.6 71.3 ± 2.5 70.2 ± 1.7 

n ED E,Bf 0.0559 ± 0.0050 0.0557 ± 0.0048 0.0537 ± 0.0032 

EDE Mocker Model 

w < -0.694 < -0.708 < -0.909 

C < 2.950 < 2.623 < 3.214 

n m 0.307 ±0.050 0.303 ±0.048 0.272 ±0.021 

a 8 0.716 ±0.072 0.714 ± 0.070 0.744 ± 0.049 

Ho 66.1 ±4.4 66.3 ±4.3 69.3 ± 1.8 

«ede(zi ss ) < 0.037 < 0.036 < 0.026 

n EDE;S f 0.0681 ±0.0144 0.0682 ±0.0139 0.0644 ± 0.0104 




0.5 0.6 0.7 0.8 0.9 55 60 65 70 75 0.04 0.06 0.08 0.1 0.12 



FIG. 17: One-dimensional marginalized likelihood constraints 
on the dark energy parameters wo and C, as well as f2 m , 
Ho, as and SlEDE.sf from different present data combinations: 
WMAP5 only (black solid lines), WMAP5±ISW (red dashed 
lines) and All datasets (blue dash-dot lines). 



more stringent. 

When combining all the datasets together, the con- 
straints tighten significantly: the error bars of some pa- 
rameters are reduced by a factor of 1.5. These is due to 
the constraining power of SNIa and BAO. These results 
are consistent with some previous works [l|, HI] and also 
with the recent findings based on the clustering of lumi- 
nous red galaxies of the DR7 Thus, present data 
sets allow to constrain the amount of dark energy in the 
structure formation era at the percent level. 



2. EDE Mocker Model 

Due to the lack of cosmological probes, the behaviour 
of the dark energy component is very poorly constrained 
in the redshift range 2 < z < 1100. Therefore, when 
using WMAP5 data only, the constraints on the param- 
eters wq and C, describing the equation of state of early 
dark energy, are very weak, namely the 95% upper limits 
are wo < —0.694 and C < 2.950. Consequently, current 
observations still allow very large amount of early dark 
energy at high redshift z ~ 1090 as ^ede^Iss) < 0.037 
(95% CLA, which is consistent with the results obtained 
by Refs. [il|, [54], [9(3]. Early dark energy models with a 
non-negligible fraction of dark energy density still fit the 
CMB data very well. 

When comparing the results of ACDM and EDE 
model, the error bars of some parameters are signifi- 
cantly enlarged by a factor of two, shown in Fig ll7l 
due to the degeneracies between the dark energy pa- 
rameters and other background parameters. We find 
that the constraint on the parameter OEDE.sf becomes 
rather weak for EDE. The current constraint on f^EDE.sf 
is 5lEDE,sf = 0.0681±0.0144 at 1 a confidence level, which 
is obviously higher than the pure ACDM model, although 
the two agree at la level: ^EDE,sf = 0.0559 ± 0.0050 
(lcr). This is because of the higher dark energy abun- 
dance at high redshifts with respect to ACDM. Moreover, 
the linear growth factor of early dark energy model is 
suppressed significantly by the large value of parameter 
C. When the fraction of dark energy density becomes 
large in the EDE model, the present value of as will 
be lower, erg = 0.716 ± 0.072 (lcr), which is obviously 
lower than one obtained in the pure ACDM framework: 
as = 0.797 ± 0.035 (1 a), while the error bar is enlarged 
by a factor of two [4l|] . In Fig. [TS] we can clearly see the 
anti-correlation between as and f^EDE,sf ■ 

Although the ISW data are directly related to the dark 
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FIG. 18: Two dimensional contour in the (ag, fiEDE.sf) panel 
from different present data combinations: WMAP5 only 
(black solid lines) , WMAP5+ISW (red dashed lines) and All 
datasets (blue dash-dot lines). 



energy parameters and contain information on the low 
redshift universe (z < 5), at present the constraints are 
too weak to offer much improvement. In fact, the re- 
sults do not improve significantly when adding the ISW 
data, which is shown in Fig. [501 and Table V. The con- 
straints on wo and C improve slightly: wq < —0.708 and 
C < 2.623 at 95% confidence level which are close to 
the pure ACDM model. Meanwhile, the constraints on 
other cosmological parameters, such as <7s and ^EDE.sf, 
also tighten a little, but not significantly. In Fig. [15] we 
plot the two dimensional constraint on (erg, ^EDE,sf ) from 
different data combinations. WMAP5+ISW data com- 
bination gives clearly tighter constraints than WMAP5 
only. 

Finally, we add some low-redshift observational data, 
such as SNIa and BAO data. Due to their constrain- 
ing power, the constraint on wq improves significantly: 
wo < —0.909 at 95% confidence level. However, the 95% 
upper limit on C has not been improved: C < 3.214, 
because of the anti-correlation between wo and C [4l| . 
The constraints on other parameters, when combining 
all datasets together, become slightly more stringent and 
are consistent with the previous work |4lj . in which a 
similar analysis was carried using gamma-ray bursts and 
Lyman-a forest observations instead of the ISW effect. 



VIII. FUTURE RESULTS 

From the results presented above, we see that both in 
ACDM and EDE model, the ISW data give a little im- 
provement on the constraints on the cosmological param- 
eter we considered, when compared to other observations, 
such as SNIa and BAO. Therefore, it is worthwhile dis- 
cussing whether future ISW data could give more strin- 



TABLE VI. Constraints on the ACDM model from the 
future measurements. Here we show the standard deviations. 



Parameter CMB CMB+ISW All Datasets 
ACDM Model 



b 






0.1943 


a 






0.0016 


■^amp 






0.2856 




0.0090 


0.0080 


0.0029 




0.0075 


0.0070 


0.0054 


H 


0.7730 


0.7012 


0.2327 


^EDE,sf 


0.0014 


0.0013 


0.0005 



gent constraints on the cosmological parameters. For 
this purpose we have performed a further analysis and 
we have chosen two fiducial models in perfect agreement 
with current data: a pure ACDM model and an EDE 
model with parameters taken to be the best-fit values 
from the current constraints of "All" datasets combina- 
tion. 



A. ACDM Model 

Firstly, we use the fiducial ACDM mock datasets to 
constrain the parameters in the ACDM model, as well 
as other three parameters, b, a and A amp . In Table 
VI we list the standard deviations of those parameters 
from these mock future measurements. We remind that 
the mock ISW data sets consist of data with covariance 
matrix reduced by a factor nine when compared to the 
present ACF and CCF data. 

Due to the smaller error bars of the mock data sets, the 
constraints on the cosmological parameters from CMB 
only improve significantly by a factor of three, when com- 
paring to the current results. When adding the simulated 
ISW data, the constraints improve further and the im- 
provements are larger than those from the current obser- 
vations, since the ISW data with smaller error bars are 
now more helpful in breaking the degeneracies between 
the parameters that keep the CMB angular diameter dis- 
tance unchanged. 

When using all datasets together, we present the con- 
straints on the cosmological parameters, as well as those 
of three parameters. The constraints on the parameters 
improve significantly. In particular, due to the more ac- 
curate ISW data, the standard deviations of b, a and 
Agjnp have been shrunk to 0.19, 0.15% and 0.29. In this 
case, the significance of ISW signal will be clearly en- 
hanced. The future CMB measurement and galaxy sur- 
vey could be very useful to detect the ISW effect at a 
much higher significance than now and also to constrain 
the bias at high redshift. 
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TABLE VII. Constraints on the early dark energy model from the future measurements. Here we show the standard 
deviations. For some parameters that are only weakly constrained we quote the 95% upper limit. 



Parameter 


CMB 


CMB+ISW 


All Datasets 


CMB 


CMB+ISW 


All Datasets 




Fiducial ACDM Model 


Fiducial EDE Mocker Model 


w 


< -0.4952 


< -0.6586 


< -0.9571 


< -0.6154 


< -0.6455 


< -0.9041 


C 


< 1.4452 


< 1.8473 


< 2.7632 


< 3.3496 


< 1.9581 


< 2.4893 




0.0642 


0.0412 


0.0031 


0.0494 


0.0406 


0.0033 




0.0607 


0.0425 


0.0108 


0.0512 


0.0433 


0.0188 


H 


5.3335 


3.7914 


0.3186 


5.0357 


3.8344 


0.4733 


^ede(zi ss ) 


< 0.0063 


< 0.0057 


< 0.0035 


< 0.0141 


< 0.0129 


< 0.0122 


^EDE.sf 


0.0055 


0.0039 


0.0018 


0.0060 


0.0054 


0.0043 
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FIG. 19: One dimensional marginalized likelihood constraints 
on the dark energy parameters wo and C, as well as f2 m , Ho, 
<Tg and f^EDE.sf from different mock futuristic data combina- 
tions with the fiducial ACDM model: CMB only (black solid 
lines), CMB+ISW (red dashed lines) and All datasets (blue 
dash-dot lines). 



the constraints on the parameters improve significantly, 
unlike for the current results. The 95% upper limit 
of current equation of state of dark energy wq is now 
w < -0.659, while w < -0.495 (95% C.L.) obtained 
from CMB data only. By contrast, the constraint on C 
has not been improved significantly: C < 1.847 (95% 
C.L.), due to the anti-correlation between wq and C. 
Meanwhile, all of the related parameters have been con- 
strained more tightly than ones from CMB data only by 
a factor of two. This simulated ISW data with smaller 
error bars is effective in breaking the degeneracies among 
the parameters. 

Furthermore, when using all datasets, from Fig.[T9"land 
Table VII we can find that the constraints of many pa- 
rameters have been tightened significantly and the de- 
generacies have been broken further. The constraint on 
wo becomes very tight w < —0.957 at 95% confidence 
level, due to the accurate Supernovae data. And since 
the 95% upper limit of Oede^iss) is also suppressed ap- 
parently, f^EDE^iss) < 0.0035, many early dark energy 
models can be ruled out. 



2. Fiducial EDE 



B. EDE Mocker Model 

1. Fiducial ACDM 

Firstly, we choose the pure ACDM as the fiducial 
model. In table VII we list the forecasts for some re- 
lated parameters using the future measurements. 

From Fig. HU and Table VII, we can find that the CMB 
data only cannot constrain the cosmological parameters 
well, as we expect. In fact, the standard deviations of pa- 
rameters are rather large. Early dark energy models with 
a non-negligible fraction of dark energy density cannot be 
ruled out by the CMB data only, namely the constraint 
on the fraction of dark energy density at z ~ 1090 is 
^ede(ziss) < 0.0063 (95% C.L.). 

Interestingly, when including the future ISW data, 



We also choose the EDE fiducial model to determine 
the cosmological parameters. The one-dimensional con- 
straints of some related parameters from different data 
combinations are plotted in Fig. [50] We obtain similar 
conclusions to the ACDM case. 

The CMB data only cannot constrain the parameters 
very well. When including the ISW data, the constraints 
on the parameters improve significantly. However, even 
using CMB+ISW data combination, the best fit values 
of wo and C are still close to the ACDM model. In this 
case, early dark energy models cannot be distinguished 
from the pure ACDM model. 

Finally, the "All" datasets combination give the most 
stringent constraints on the parameters. In this case, the 
peaks of one dimensional distributions of wq and C are 
moving away from the ACDM model, too = — 1, C = 0, 
see Fig. [20] The 95% confidence level are w < -0.904 
and C < 2.489, respectively. These results imply that 
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FIG. 20: One dimensional constraints on the dark energy 
parameters wo and C, as well as f2 m , Ho, as and SlEDE,sf from 
different mock futuristic data combinations with the fiducial 
EDE model: CMB only (black solid lines), CMB+ISW (red 
dashed lines) and All datasets (blue dash-dot lines). 



the future measurements could distinguish between the 
pure ACDM model and early dark energy models. 



IX. CONCLUSIONS AND DISCUSSIONS 

Most cosmological probes are sensitive to the be- 
haviour of dark energy at very low redshift; thus, it is 
of great interest to investigate those observable which 
are able to complement the analysis with constraints on 
the dark energy abundance at high redshifts, close to the 
onset of cosmic acceleration and possibly to the regime 
which is probed by other large scale structure observ- 



ables like the Lyman-a forest [9l|. In this paper we 
exploited the capabilities of the ISW effect in the high 
redshift universe, using the cross-correlation signal be- 
tween the one million photometrically selected QSOs of 
the SDSS DR6 catalog El and the CMB maps of the 
WMAP year 5 satellite (H!|. From the SDSS DR6-QSO 
catalog we extract the QSO ACF and estimate the bias 
and stellar contamination by considering several redshift 
ranges between z = and z = 5. We have given par- 
ticular emphasis to our subsamples at z > 1.5 where an 
overall weak evidence for a non-zero ISW signal is found 
at the 1.5cr confidence level. The evidence is instead 
at the 2.7 a level if we consider the whole sample and 
this is in agreement with investigations based on similar 
samples |14l Il5j . Our high-redshift sample has a mean 
redshift of z — 2 and z — 2.3 if we choose the conserva- 
tive QSO sample or the speculative one (sources with no 
ultra-violet excess), respectively. This is a new regime 
when compared to other probes that have been more ex- 
tensively used for ISW studies such as galaxies or X-ray 



observations. 

We have focused on modifications to the standard 
ACDM cosmology that results in a non-negligible amount 
of dark energy in the structure formation era: these mod- 
els are generically labeled as Early Dark Energy models 
and are characterized by an early departure for matter 
dominance in the cosmic expansion, with consequences 
for background evolution and therefore structure forma- 
tion. We found that at present the results are still rather 
weak to provide competitive constraints on the parame- 
ters describing either the ACDM or EDE models, even 
though constraints at the percent of sub-percent level on 
the energy density contribution of this component in the 
structure formation era can be achieved. Adding present 
high redshift ISW data to CMB data from WMAP does 
not improve the constraints significantly. However, we 
also forecast future performance of QSO data by assum- 
ing a measurements of CCF and ACF data by reducing 
the error bars by a factor three and combining these with 
mock Planck data. In this case the results are partic- 
ularly interesting since the improvement to CMB data 
alone when adding the ISW information can result in a 
factor 1.5 on most of the cosmological parameters. If we 
further add some Supernovae luminosity distance moduli 
like those that can be provided by the SNAP satellite the 
constraints can become even tighter and up to a factor 
between 3-10 better than the CMB alone. 

Here we summarize our main conclusions in more de- 
tail: 

• We compute the QSO Auto-Correlation Function 
(ACF) from SDSS DR6-QSO and extract bias and 
catalog efficiency: we measure the quasar-matter 
bias and the typical error on the bias is at the level 
of Ob = 0.8 when we split in five bins the sam- 
ple and of Ob = 0.3 when we consider the whole 
sample. We present results for a more conserva- 
tive selection of QSOs using also ultra-violet excess 
flags and a less conservative selection that consider 
all the sources of the catalog. 

• EDE models usually result in higher values for the 
bias than those of ACDM by ~ 10%, due to the 
slower growth of density perturbations in EDE cos- 
mology that requires a higher bias to match the 
observed value. 

• We compute CCF values by cross-correlating the 
QSO number density with the CMB temperature 
and found an evidence at 2.7 a level for the all sam- 
ple for an ISW effect. At high redshift this evidence 
reduces to 1.5 a, which interestingly compares with 
the prediction of null detection in ACDM. This 
non-vanishing signal could be caused by the large 
bias, stellar contamination and large number of 
quasar sample, or an higher dark energy abundance 
at high redshifts. However, at present this signal 
is still too weak to distinguish between the ACDM 
and early dark energy models. 
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• The parameter A amp was used to quantify the 
disagreement between theoretical predictions and 
data. We found that this disagreement vanishes 
in the high redshift subsample, either for EDE or 
ACDM models and is instead at the ~ 2 a level for 
the low-redshift subsamples. 

• We forecasted futuristic large-scale structure mea- 
surement of the ISW effect using the cross- 
correlation that can be provided by photometric 
surveys. Assuming smaller error bars by a factor 
of three, the more accurate ISW data is helpful to 
break the degeneracies among some parameters and 
results in an overall improvement on mock Planck 
data alone by a factor 1.5. Furthermore, the stan- 
dard derivations of those three parameters, b, a 
and A amp , arc also shrunk significantly. The fu- 
ture CMB measurement and galaxy survey could 
be very useful to detect the ISW effect at a much 
higher significance. 

• Further adding Supernovae measurements improve 
results by a factor 3-10 due to the better determina- 
tion of H and f2 m that can be provided by SNAP 
observations. 

The use of this state-of-the-art QSO catalogue can 
open up a completely new window on the high-redshift 
large scale structure of the universe and allowing for a 



quantitative use of the ISW effect in a regime z > 1.5, 
which is at present weakly probed by observations. Ei- 
ther a confirmation of the ACDM model or possible de- 
partures induced by modified gravity or quintessence 
could be of fundamental importance and will be ad- 
dressed by the surveys of the near future. 
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